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Searches for Mixing and CP Violation in the D^-D^ System 

A. J. Schwartz 

Department of Physics, University of Cincinnati, P.O. Box 21 0011, Cincinnati, Ohio 45221, USA 



We review searches for mixing and CP violation in the D^-D " system and discuss the first evidence for mixing 
recently obtained by the Belle and BABAR collaborations. We also present world average values for the mixing 
parameters x = Am/F and y = Ar/(2r) as calculated by the Heavy Flavors Averaging Group. 
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1. Introduction 

Mixing in the D^-D " system has been searched for 
for more than two decades without success — until 
this past year. Both "B-factory" experiments, BeUe 
and BABAR, have recently published evidence for this 
phenomenon [1, 2]. Here we review these measure- 
ments and discuss their implications. In total four 
measurements are presented, involving the following 
decay modes [3]: D^^K+l'v, D" K+R- /n+n- , 
D° -^K+n-, and D° tt+tt- . 

Mixing in heavy flavor systems such as that of 
and B^ is governed by the short-distance box dia- 
gram. However, in the system this diagram is both 
doubly-Cabibbo-suppressed and GIM-suppressed rel- 
ative to the amplitude dominating the decay width, 
and thus the short-distance rate is very small. Conse- 
quently, D^-IJO mixing is expected to be dominated 
by long-distance processes that are difficult to calcu- 
late; theoretical estimates for the mixing parameters 
X = Am/r and y = Ar/(2r) range over two-three 
orders of magnitude [4]. Here, AAf and AF are the 
differences in the masses and decay widths, respec- 
tively, of the two D^-D^ mass eigenstates, and F is 
the mean decay width. 

All methods discussed here identify the flavor of the 
/D^ when produced by reconstructing the decay 



D* 



D°TT+ or D* 



^ D 't:^; the charge of the 
accompanying pion identifies the D flavor. Because 
M^. -Mjjo-M^^ « 6 MeV, which is relatively small, 
the pion has very low momentum and thus is denoted 
TTg ("tt slow"). The decay time {t) is calculated 
via {i/p) X Mjj, where £ is the distance between the 
D* and decay vertices and p is the D'^ momentum. 
The D* vertex position is taken to be the intersection 
of the TTj, momentum with the beamspot profile. Most 
of the precision on t is due to the vertical {y) compo- 
nent of £, as the spread of the beamspot is only a few 
microns in this dimension. To reject D^*) decays orig- 
inating from B decays, one requires p^* > 2.5 GeV, 
which is the kinematic endpoint. 



2. "Wrong-sign" D^^K+l'i^ Decays 

The decay K^i^v is in principle an ideal 

signature for mixing, as this "wrong-sign" (WS) fi- 



nal state can be reached only via a D°^D° transi- 
tion (in contrast to the "right-sign" (RS) final state 
K^i^u). However, the neutrino in the final state pre- 
cludes the decay from being fully-reconstructed, and 
consequently the decay time measurement is smeared. 
The decay time dependence is given by e~'"*[l + {x^ + 
y^)t^/2], and both Belle and BABAR measure the co- 
efficient (x^ + y'^)/2 = Rj^ by fitting the decay time 
distribution. To reduce backgrounds, only the elec- 
tron channel [Kev) has been used. 

Belle uses 253 fb~^ of data [5], making a relatively 
loose event selection and fitting the resulting M{Kev) 
and AM = M{-k^ Kev)-M{Kev) distributions to de- 
termine the yield of WS events. To improve the resolu- 
tion on both AM and the decay time, two corrections 
are made to better determine P^, = -Pcms^^Trife^^rcst j 
where P^est is the total four-momentum of all tracks 
and photons in the event besides the tt^, K, and e 
candidates, and P^^^ is the four-momentum of the 
center-of-mass (cms) system. First, P^.^^^ is rescaled 
by a factor x such that (P^^^-xP^^^^)"^ — P|,. = M|,. ; 
second, p^^st is rotated in the plane of Pj.^^^ and p^ 
such that (P„„, - P^^^^ - xP„,t)2 ^ p2 ^ 0. 

The ratio of WS to RS events are measured for six 
bins of decay time ranging from one to ten life- 
times. Summing over all bins, 16430 ± 155 RS and 
— 1 ± 37 WS events are found. For each bin the ratio 
depends on P^^, and the six resulting values of Rj,^ 
are fit to a constant. The result is 



R 



(0.020 ±0.047 ±0.014)%, 



(1) 



where the first error listed is statistical and the second 
error is systematic. The latter is dominated by uncer- 
tainty in the background A A/ PDF. As much of the 
P^ likelihood is in the unphysical (negative) region. 
Belle uses a Feldman- Cousins [6] approach to calcu- 
late an upper limit; the result is Rj^.j < 0.10% at 90% 
C.L. This limit implies both |x| and \y\ are < 4.5%. 

The BABAR experiment uses 344 fb"^ of data [7] 
and, in contrast to Belle, imposes tight selection crite- 
ria to reduce background as much as possible (with a 
corresponding loss in efficiency) . The most restrictive 
criterion is that events must have a Z)'*) decay fully 
reconstructed in the hemisphere opposite that of the 
semileptonic decay, where the modes D'^tt'^, 
D° K-n+, D" K-n+TT°, D° K-tt+tt+tt-, 
and Z)+ ^ if ^7r+7r+ arc used. This "double-tagging" 
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eliminates copious background from real D^'s combin- 
ing with random tt^ tracks to make false D* candi- 
dates, but it reduces the signal efficiency by an order 
of magnitude. The neutrino momentum is determined 
via a neural network algorithm, and events are re- 
quired to have decay times in the range 600-3900 fs 
(corresponding to 1.5-9.5 -D" lifetimes). A separate 
neural network is used to select signal events, and a 
final set of kinematic selection criteria are applied to 
the signal side. These criteria include cuts on the ir^ 
and e momenta, and on the dE/dx of the electron 
track in the silicon vertex tracker. 

The final candidate samples for WS data and Monte 
Carlo (MC) simulation are shown in Figs, la and lb, 
respectively. In the data, three events arc observed 
in the signal region AM < 0.20 GeV/c^, whereas the 
MC predicts (for the luminosity of the data) 2.85±1.43 
background events. Together these values give = 
(0.004 

^aoe)^! where the error corresponds to where 
the log-likelihood function for (the true number of 

signal events) rises by 0.50 units with respect to the 
minimum value. The points where the log-likelihood 
function rises by 1.35 units give a 90% C.L. constraint 
-0.13% < Rj^ < 0.12%. This upper bound is similar 
to the upper limit obtained by Belle. 
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Figure 1: WS AM distribution for BABAR's D° 

candidate sample [7], in data (left) and MC (right). 
The MC sample is 1.75 times the size of the data sample. 
The dark (light) histogram shows events after (before) the 
final kinematic selection. In the left plot, region "1" de- 
notes the signal region used to determine R^- 



3. Hadronic Decays -D° h+ h'^ 

Two-body hadronic final states K~^K~ 
K^TT^ can be reached from either or D°; thus 
D^^f and ^ f amplitudes both contribute 

to the decay rate, and detecting the effect of the latter 
provides evidence for mixing. The time dependence of 
the decay rate TZ is given by 



-rt 



(y ReA — xlmX)t ■ 



where A = {q/p)A{D" I) / A{D° /) and q,p 
are complex coefficients relating mass eigenstates to 



flavor eigenstates: 2 = P\D^) ± q\D^). The pa- 
rameter A can be written —\q/p\\Aj^^J Ajja\e'^'''^^^\ 
where 5 is the strong phase difference between ampli- 
tudes A{D^^f) and A{D°^f), and ^ is a possible 
weak phase. In the absence of CP violation {CPV), 
\q/p\ = 1 and (/) = 0. For / = K-tt+, A(D° ^ f) is 
doubly-Cabibbo-suppressed, A{D^ f) is Cabibbo- 
favored, and thus |A| -C 1 and TZ » e~^*. 



3.1. CP-eigenstates K+K~ and tt+tt" 

For decays to self-conjugate states K'^K~ and 
TT+TT", |A| = \q/p\ ~ 1 and the third term in Eq. (2) 
can be neglected since and \y\ are very small. 
As 6 — 0, if there is no CPV in mixing {\q/p\ = 1) 
then A = -e"^ and TZ 

„-r(l+ycos0)t 



observable r 



K 



"h+h 



Thus t^-^+/t; 



e-"[l - {y cos (t>)Tt\ » 

~1 + ?/cos(/). The 
1 is denoted y^p and, for 



no CPV in mixing, equals y cos 0; if CP is conserved, 
Vcp — y- Allowing for arbitrary CPV, one obtains [8] 



ycp = \ 



2/ cost 



a;sin( 



The (normalized) difference 
lated expression [8] 



in lifetimes Ay 



is equal to the re- 



y cos (p — 



+ 



a;sin( 



This method has been used by numerous exper- 
iments to constrain y^jp [9]. Belle's new measure- 
ment [1] uses 540 fb~^ of data and both K'^K~ and 
TT+TT" final states. One advance of this analysis is the 
resolution function, which is constructed as a sum over 
many Gaussian functions G: 



k=l 



with standard deviations a- 



pull 



X cr. 



this expression, is the weight of the value taken 
from the normalized, binned, — > K~n~^ distribu- 
tion of (Tf, the event-by-event uncertainty in the de- 
cay time (i.e., = crj(bin i)). Parameter is the 
weight of value cr^"^' obtained by fitting the MC pull 
distribution to a sum of three Gaussians with widths 
cr^"" (fc = 1 — 3) . The are scale factors to account 
for differences between MC and data, and is a com- 
mon offset. The parameters Sj. and are left free 
when fitting for ^.^p. This resolution function, and 
a slight variation with an additional offset parameter, 
) (2) yields accurate values of the D° ^K~tt~^ lifetime over 
all running periods. The mean value is 408.7 ± 0.6 fs, 
which is consistent with the PDG value [10] (and ac- 
tually has greater statistical precision). 
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Fitting the K^tt^, K^K~, and tt+tt" decay time 
distributions (Figs. 2a-c) shows a statistically signif- 
icant difference between the K~'jt~^ and life- 
times. The effect is visible in Fig. 2d, which plots 
the ratio of event yields iV^+^./A'^^^ as a function 
of decay time. Performing a simultaneous maximum 
likelihood (ML) fit to all three hfh2 samples gives 



Vcp = (1.31 ±0.32 ±0.25)%, 



(3) 



which deviates from zero by 3.3o-. The systematic 

error is dominated by uncertainty in the background 
decay time distribution, variation of selection criteria, 
and the assumption that is equal for all three final 
states. The analysis also measures 



(0.01 ±0.30 ±0.15)%, 



(4) 



which is consistent with zero (no CPV). The sources 
of systematic error for Aj. are similar to those for y^p. 




Figure 2: Projections of the decay-time fit superimposed 
on the data for D° K+K~, D° K-n+, and D° 
tt+tt" decays, from Belle [1]. The hatched area represents 
the background contribution. Plot (d) shows the ratio of 
decay-time distributions for (K'^ K~ + ■k'^'k~) and 

D'^ -^K~tv'^; the solid line is a fit to the points. 



3.2. "Wrong-sign" D^^K+tt- Decays 



For / = if 
fies to e-r* [R 
R 



D 




CPV, Eq. (2) simpli- 
{x'^ + y'2) t'^/4\ , where 



)/A{D° ^ K+TT-) and x' = 
xcos6 + y Hill 6. y' = ycoaS — a; sin 5 are "rotated" 
mixing parameters. Both BABAR [2] and Belle [11] 
do unbinned ML fits to the decay-time distribution of 



WS D^^if+TT decavs to determine 



! y 



and Rr 



small and the mixing terms in the above expression 
play a larger role; however, there is substantial back- 
ground, ~48%. The largest background component 
consists of real D'^ K^tt^ decays combining with 
random tt" tracks; fortunately, the decay time distri- 
bution for this background is simple, the same as that 
for RS D°^K-Tr+ decays. 

The results of the BABAR and Belle fits are listed in 
Table L For BABAR, the value for x''^ is negative, i.e., 
outside the physical region, nominally due to statisti- 
cal fluctuation. The BABAR likelihood contours are 
shown in Fig. 3. The no- mixing point (0, 0) has AC = 
23.9 units above the minimum value, corresponding 
to a CL of only 0.01% (3.9it) including systematic un- 
certainty. This constitutes evidence for mixing. The 
largest systematic error is from uncertainty in mod- 
eling the tail of the background decay time distri- 
bution. The mixing is visible in Fig. 4, which plots 
the ratio of the background-subtracted yields of WS 
to RS decays in bins of decay time. For each bin, 
the yields are determined from two-dimensional fits 
to variables M{Ktt) and AM = M(7r^ K'k)-M{Ktt). 
The plot shows the ratio increasing with decay time, 
consistent with Eq. (2) but inconsistent with the no- 
mixing or flat hypothesis. Fitting to Eq. (2) gives 
X^/dof = 1.5, whereas fitting to a flat distribution 
gives xVdof = 24.0. To allow for CPV, BABAR fits 
the D° and D° samples separately; the results are 
consistent with each other, showing no evidence of 
CPV (see Table I). 



Table I BABAR [2] and Belle [11] results from fitting the 

decay time distribution of D*^ —>K^tt~ decays. The errors 
listed are statistical plus systematic, except for those from 



BABAR's D" 
onh'. 



and D subsamples, which are statistical 



Exp. (fb^^) 


x''' (%) 


y' (%) 


Rd (%) 


BABAR (384) 
D° only 
D° only 


-0.022 ± 0.037 
-0.024 ± 0.043 
-0.020 ± 0.041 


0.97 ± 0.54 
0.98 ± 0.64 

0.96 ±0.61 


0.303 ±0.019 


Belle (400) 
CPV-allwd 


0.018 1'>:;- 

< 0.072 


06 

^•'J" -0.39 

(-2.8,2.1) 


0.364 ± 0.017 



Because JsT+tt is doubly-Cabibbo-suppressed, Rj^ is 



The Belle measurement has somewhat greater sta- 
tistical precision than that of BABAR, but the central 
values arc in the physical region (x''^ > 0). Belle ob- 
tains confidence regions for x''^ and y' using a toy-MC 
frequentist method. Due to the proximity of the un- 
physical region, the procedure uses Feldman-Cousins 
likelihood ratio ordering [6]. The resulting contours 
arc shown in Fig. 5. The CL of the no- mixing point 
(0,0) is 3.9%, corresponding to 2.1cr. The largest sys- 
tematic uncertainty arises from variation of the Pjj, 
minimum value cut. 

Belle searches for CPV by fitting the D*^ and 
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Figure 3: Two-dimensional la-5a contours for (x'^,2/'), 
from BABAR^s fit to the decay time distribution of D° — > 
K^-K~ decays [2]. The contours are calculated from the 
change in the value of — 21n£ from the minimum value, 
and systematic uncertainties are included. The point is 
the best-fit value. 
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unmeasured). The fit results are (f) = (9.4 ± 25.3)° or 
(84.5 ± 25.3)° for the same or opposite signs of 
and x'-, and Aj^ = 0.67 ± 1.20. 
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Figure 5: 95% C.L. contours for {x'^, y') from Belle's anal- 
ysis of — > it'+TT" decays [11]: dotted (dashed) is statisti- 
cal (statistical plus systematic) contour for no CPV; solid 
is statistical plus systematic contour allowing for CPV . 
The point is the best-fit value for no CPV . The contours 
are obtained from a frequentist toy-MC calculation. 



4. Dalitz Plot Analysis of 



+7r- 



Figure 4: The WS K'*'n~ branching fraction for 

separate bins of decay time, from BABAR [2]. The dashed 
line shows the expectation based on the values of {x'^,y') 
obtained from the decay time fit. The data agrees well 
with the dashed line but not with a fiat (i.e., no mixing) 
line. 



The time dependence of the Dalitz plot for D''^ 

decays is sensitive to mixing parameters x 
and y without ambiguity due to strong phases. For 
a particular point in the Dalitz plot {rri^, m?_), where 
m_|_ = m{Kg'K'^) and m_ = m{K^iT~), the overall 
decay amplitude is 



distributions separately. 
{x'^+,y'+,Ro) and [x'' 



The results, denoted 
i?^), respectively, are 
used to calculate the CPV parameter Aj^ = {R^ — 

The- 



Rm)I{Rm + Rm)^ where R 



'■M 



(a;'±2 + y'±2)/2. 



oretically, Aj^ = {\q/p\^ -l)/{\q/p\'^ + l) « \q/p\''-l. 
The CPV parameter ATg{q/p) = <p and mixing pa- 
rameters x'^, y' are determined via the relations 



X- = [{l±Aj^)/{lT A^)Y'\x' cos cj)±y' sin, 
[(1 ± Am) /{I T Aj^)]"\y' cos ^Tx' sine 



The resulting confidence region for {x'^,y') is plotted 
in Fig. 5 as the solid contour; the complicated shape 

is due to there being two solutions for {x\y'), de- 
pending on the relative sign of x'^ and x'~ (which is 



Ajjo (m_|_,m 



2 2s('liJ)+<'2{t) 



2/2)« 



The first term 



where e^^ 2){t) = e 
represents the (time-dependent) amplitude for — > 
ii'^ TT+TT", and the second term represents the ampli- 
tude for £)0^Z)"^ifg7r+7r-. Taking the modulus 
squared of Eq. (5) gives the decay rate or, equiva- 
lently, the density of points p{rn'^, m,'^;t). The result 
contains terms proportional to cosh(yrt), cos{xrt), 
and sin(a;rt), and thus fitting the time-dependence of 
p{m\,m?_;t) determines x and y. This method was 
developed by CLEO [12]. 
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To use Eq. (5) requires choosing a model for the 
decay amplitudes -4^^ — ^(m^, m^). This is usually 

taken to be the "isobar model" [13], and thus, in ad- 
dition to X and y, one also fits for the magnitudes 
and phases of various intermediate states. Specifi- 
cally, Ajjo irn\ , mt ) = J2j ^^^^ > where Sj is a 
strong phase, Aj is the product of a relativistic Breit- 
Wigner function and Blatt-Weiskopf form factors, and 
the parameter j runs over all intermediate states. This 
sum includes possible scalar resonances and, typically, 
a constant non-resonant term. For no direct CPV, 
A—g{m'^,m^) = .4^0 m^); otherwise, one must 

consider separate decay parameters (a^ , S^) for de- 
cays and {(ijjSj) for decays. 

Belle has recently fit a large D'^ tt^7t~ sample 

selected from 540 fb""'^ of data [14]. The analysis pro- 
ceeds in two steps. First, signal and background yields 
are determined from a two-dimensional fit to vari- 
ables MiKiTTT) and AM = M{iT^KinT) - M{KinT). 
Within a signal region \M{KTnT) - Mj^o \ < 15 MeV/c^ 
and \AM - 5.9 MeV] < 1.0 MeV (corresponding to 
3(7 in resolution), there are 534000 signal candidates 
with 95% purity. These events are fit for x and y; 
the (unbinned ML) fit variables are m^, 77i?_, and 
the decay time t. Most of the background is com- 
binatoric, i.e., the candidate results from a ran- 
dom combination of tracks. The decay-time distribu- 
tion of this background is modeled as the sum of a 
delta function and an exponential function convolved 
with a Gaussian resolution function, and all parame- 
ters are determined from fitting events in the sideband 
30 MeV/c^ < iMiKiTTT) - M^oj < 55 MeV/c^. 

The results from two separate fits are listed in Ta- 
ble II. In the first fit CP conservation is assumed, i.e., 
q/p=l and A—„{m'^,mt) = Ajjoi'm'^,m\). The free 
parameters are x, y, t^o , some timing resolution func- 
tion parameters, and decay model parameters (a^, S^.). 
The results for the latter are listed in Table III. The 
results for x and y indicate that x is positive, about 2a 
from zero. Projections of the fit are shown in Fig. 6. 
The fit also yields r^, = (409.9 ± 0.9) fs, which is 
consistent with the PDG value [10] (and actually has 
greater statistical precision). 

For the second fit, CPV is allowed and the 
and samples are considered separately. This in- 
troduces additional parameters \q/p\, Arg{q/p) = cj), 
and (cipSj). The fit gives two equivalent solutions, 
{x,y,(j)) and (— x, — y, -I- tt). Aside from this pos- 
sible sign change, the effect upon x and y is small, 
and the results for \q/p\ and are consistent with 
no CPV. The sets of Dalitz parameters {a^,5j.) and 
(a^jS^) are consistent with each other, indicating no 
direct CPV. Taking aj = aj and ^j=Sj (i.e., no direct 
CPV) and repeating the fit gives \q/p\ = 0.95 to;2o 
and0= (-21^?)°. 

The dominant systematic errors are from the time 



Table II Fit results and 95% C.L. intervals for x and y, 
from Belle's analysis of _D" iCg tt+tt" decays [14]. The 
errors are statistical, experimental systematic, and decay- 
model systematic, respectively. 



Fit Param. 




Result 


95% C.L. inter. 


No X (%) 


0.80 


1 n 9Q +0.09+0.10 
^ U.ztJ _Q 07 -0.14 


(0.0,1.6) 


CPV y (%) 


0.33 


1 n 9/1 +0.08+0.06 
^ '-'•^^ -0.12 -0.08 


(-0.34,0.96) 


CPV X (%) 


0.81 


1 n +0.10+0.09 

^ U.OU _o.07 -0.16 


\x\ < 1.6 


y(%) 


37 +n 25+°-°''+'' °'' 

U.OI ^ U.ZO _Q _Q Qg 


\y\ < 1.04 


\<1/P\ 


0.86 


+0.30+0.06 1 n n« 

-0.29 -0.03 ^ 




4>n 




1/1 +16+5+2 
-18 -3 -4 






ml (GeV^/O 



(GeV^/c") 




2000 4000 
Proper time (fs) 



Figure 6: Projection of the unbinned ML fit superimposed 
on the data for ^ tt^tt" decays, from Belle [14]. In 
(d), the hatched area represents the combinatorial back- 
ground contribution, and the lower plot shows the ratio of 
decay-time distributions for events in the i('*(892)"'" and 
K*{dt92)~ regions, where sensitivity to {x,y) is highest. 



dependence of the Dalitz plot background, and the 
effect of the pjj, momentum cut used to reject D*^s 
originating from B decays. The default fit includes tttt 
scalar resonances Ui and (T2 ; when evaluating system- 
atic errors, the fit is repeated without any tttt scalar 
resonances using JsT-matrix formalism [15]. The infiu- 
ence upon x and y is small and included as a system- 
atic error. 

The 95% C.L. contour for [x, y) is plotted in Fig. 7. 
The contour is obtained from the locus of points where 
— 21n>C rises by 5.99 units from the minimum value; 
the distance of the points from the origin is subse- 
quently rescaled to include systematic uncertainty. 
We note that for the CPF-allowed case, the reflec- 
tions of the contours through the origin are also al- 
lowed regions. 
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Table III Fit results for D" -* K^tt+tj-- Dalitz 
plot parameters, from Belle [14]. The errors are 
statistical only. The fit fraction is defined as the 
ratio of the integral J [a^ Ar{rri^ ,m\)\^ drn^ dm\ to 



Resonance 


Amplitude 


Phase (deg) 


Fit fraction 


K*{892)- 


1.629 ± 0.006 


134.3 ±0.3 


0.6227 


7^0 (1430) - 


2.12 ± 0.02 


-0.9 ±0.8 


0.0724 


7i-2*(1430)- 


0.87 ±0.02 


-47.3 ± 1.2 


0.0133 


^"(1410)" 


0.65 ± 0.03 


111 ±4 


0.0048 


K*{1680)- 


0.60 ±0.25 


147 ± 29 


0.0002 


K*{892)+ 


0.152 ±0.003 


-37.5 ± 1.3 


0.0054 




0.541 ±0.019 


91.8 ±2.1 


0.0047 


Js:|(1430)+ 


0.276 ±0.013 


-106 ±3 


0.0013 


js:*(i4io)+ 


0.33 ± 0.02 


-102 ±4 


0.0013 


Ji'*(1680)+ 


0.73 ±0.16 


103 ± 11 


0.0004 


p(770) 


1 (fixed) 


(fixed) 


0.2111 


w(782) 


0.0380 ± 0.0007 


115.1 ± 1.1 


0.0063 


/o(980) 


0.380 ± 0.004 


-147.1 ± 1.1 


0.0452 


/o(1370) 


1.46 ±0.05 


98.6 ± 1.8 


0.0162 


/2(1270) 


1.43 ± 0.02 


-13.6 ± 1.2 


0.0180 


p(1450) 


0.72 ± 0.04 


41 ± 7 


0.0024 


Ul 


1.39 ±0.02 
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Figure 7: 95% C.L. contours for (x, y) from Belle [14]: dot- 
ted (solid) is statistical (statistical plus systematic) con- 
tour for no CPV\ dashed-dotted (dashed) is statistical 
(statistical plus systematic) contour allowing for CPV . 
The point is the best-fit value for no CPV . 
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Figure 8: World average value for j/cpi from HFAG [16]. 



5. Combining all Measurements 

All mixing measurements can be combined to ob- 
tain world average (WA) values for x and y. The 
Heavy Flavor Averaging Group (HFAG) has done 
such a combination by adding together log-likelihood 
functions obtained from analyses of 1)° ~> K^lv^ 
D" ^ h^h-, £>" ^ X+TT-, ZJO ^ X+tt-ttO, £)" -> 
X+TT^TT+TT", and 1?° — > 7r+7r~ decays, as well as 
CLEOc results for double-tagged branching fractions 
measured at the ■0(3770) resonance [16]. The com- 
bination of likelihood functions preserves correlations 
among parameters and also accounts for non-Gaussian 
errors. When using this method, HFAG assumes neg- 
hgible CPV. 

As a first step, WA values for y^p and R^^j are cal- 
culated by taking weighted averages of independent 
experimental measurements - see Figs. 8 and 9. These 
results are then converted to three-dimensional likeli- 
hood functions for [x, y, S). For example, the measure- 
ment of y(jp gives a parabolic log-likelihood function 
in y and flat distributions in x and S. The Rj^ like- 
lihood function is an annulus in the x-y plane and a 
flat distribution in 6. 

The logarithm of the likelihood functions are added, 
and the result is added to the {x, y, S) log-likelihood 
function obtained from 1)° K^tt^ decays. The lat- 
ter is determined as follows. The experiments directly 
measure a likelihood function C{x''^ , y' , Rjj)', thus one 
first projects out i?^ by allowing it to take, for any 
(x'^, y') point, its preferred value. The resulting likeli- 
hood for (x'^, y') is converted to C{x, y, 5) by scanning 
values of {x,y,5), calculating the corresponding val- 
ues of (x'^,y'), and assigning the likelihood for that 
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Figure 9: World average value for Rjij, from HFAG [16] 
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Figure 10: Two-dimensional la-5a contours for {x,y), 
obtained by adding log-likelihoods from measurements 



of D 







K+lv, D° ^ h+h-, D° 



K+-K-, D" 



K+n-n°, D° -> K+n' n+n-, and D° ifg tt+tt' de- 
cays, and double-tagged branching fractions measured at 
the VCSTTO) resonance (from HFAG [16]). 



{x''^,y') bin. This method ignores unphysical (nega- 
tive) values of x''^. The resulting function C{x,y,S) 
is added to those obtained from y^p, Rmj ^^'^ other 
measurements. The final likelihood function is pro- 
jected onto the (x, y) plane by letting 5 take, for any 
{x,y) point, its preferred value. This projection is 
shown in Fig. 10. The unusual shape around x — y — 
is mainly due to D'~' A'+tt^ decays, which disfavor 
the no-mixing point. At x = y = 0, —2\iiC rises by 
37 units above the minimum value; this difference im- 
plies that the no-mixing point is excluded at the level 
of 5.7a. 

The likelihood function is condensed to one dimen- 
sion by letting, for any value of x (y), the parameter y 
(x) take its preferred value. The resulting likelihoods 
for (x, y) give central values and 68.3% C.L. intervals 

X = (0.87;oj^)% (6) 
y = (0.66i°:^J)%. (7) 

The former is 2.6(t from zero, and the latter is 3.2(T 
from zero. 

In summary, we conclude the following: 

• the experimental data consistently indicates 
that D^'s undergo mixing. The effect is presum- 
ably dominated by long-distance processes, and 
unless \x\ 3> jyj, it may be difficult to identify 
new physics from mixing alone. 

• Since j/^p is positive, the CP-even state is 
shorter-lived, as in the K^-K^ system. How- 
ever, since x appears to be positive, the CP-even 
state is heavier, unlike in the K'^-K'^ system. 



• There is no evidence yet for CPV in the D -D 
system. 
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